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Abstract. 

Quantum systems exist at finite temperatures and are likely to be disordered 
to some level. Since applications of quantum information often rely on 
entanglement, we require methods which allow entanglement measures to be 
calculated in the presence of disorder at non-zero temperatures. We demonstrate 
how the disorder averaged concurrence can be calculated using thermal many- 
body perturbation theory. Our technique can also be applied to other 
entanglement measures. To illustrate, we find the disorder averaged concurrence 
of an XX spin chain. We find that concurrence can be increased by disorder in 
some parameter regimes. 
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1. Introduction 

Disorder is an unavoidable feature of many-body systems [Tj[2j[3]. Since the properties 
of large disordered systems are difficult to study, tools have been developed to tackle 
them such as averaging over the disorder using sampling or perturbation theory, 
or using renormalisation group techniques. For example, strong disorder in XXZ 
spin chains can be investigated using renormalisation groups [H [5]- Entanglement 
is another important feature of many-body systems [B] , one that has been shown to 
be a useful resource in many quantum information and computation schemes. It is 
therefore important to consider how entanglement in real, finite temperature systems 
is affected by disorder. 

In the context of entanglement, average disorder in spin chains has been studied 
previously El [101 HU Q21 H3 E]. This research concentrated on using sampling, 
or on renormalisation groups at zero temperature. In |15[ 116] . a perturbative 
technique to calculate a disorder averaged (finite temperature) entanglement witness 
was introduced. In this paper, the perturbative method is again used, and it is 
demonstrated how entanglement measures, rather than a witness, can be calculated. 
While both measures and witnesses are useful, considering entanglement measures 
allows any changes in the amount of entanglement by the disorder to be found. The 
disorder averaged concurrence along with other entanglement measures for weakly 
random systems are considered. 

A physical realisation of a spin chain is highly unlikely to be free from random 
variations of couplings or fields. Due to difficulties in cooling a system to its ground 
state, it is also likely to be in a thermal, mixed state. Spin chains have been shown to 
be good candidates for quantum wires, and thus have been studied extensively [T7 l fT8 ] . 
The effects of a finite temperature 1!). 20. 21 have been considered, as have the effects 
of disorder at zero temperature (HI [231 HH HSJ [2H1 [21] ; however the combined effects 
of disorder and finite temperature have not. This paper proposes how entanglement 
measures can be calculated for a random thermal spin chain. 

Two distinct averages can be taken over the disorder: quenched or annealed [3j. 
Each is useful in differing situations, though the quenched average is often considered 
the appropriate one [27) . The quenched average corresponds to the disorder effectively 
being time independent; the disorder is not in equilibrium with the system, and while 
the system evolves, the disorder remains frozen. Calculation- wise, this requires the 
thermal average to be taken before the average over the disorder, and thus we must 
average over the logarithm of the partition function. In contrast, the annealed average 
requires the average to be taken over the partition function itself. In this case, the 
disorder changes quickly, and is in equilibrium with the system as it evolves. Thus 
we take the thermal average and disorder average at the same time. Of course, the 
quenched and annealed averages are the two extremes; disorder with aspects of both 
can exist. 

Many-body perturbation theory allows us to calculate disorder averaged 
correlation functions which are crucial for quantifying entanglement. As a consequence 
of the Jordan- Wigner transformation which transforms qubits into spinless fermions, 
we can use fermionic perturbation theory to study such spin systems |28j . In this 
paper, we calculate thermal disorder averaged correlation functions and use them to 
find the disorder averaged concurrence. 

Since the quenched or annealed average over disorder can be characterised by 
taking, respectively, an average over the logarithm of the partition function or over the 
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partition function itself, we use the partition function as a starting point to calculate 
the disorder averaged concurrence. Casting the partition function into functional 
(or path) integral form, introducing a generating functional term, and replicating it, 
we find we can take the average over the disorder, and later calculate a perturbative 
expansion of the disorder averaged correlation functions for both types of average [lj[2]. 
Thus rather than calculating a direct average over the concurrence, we construct it 
from disorder averaged correlation functions. 

In particular, we consider an XX spin chain with a random term in the 
thermodynamic limit. In addition to the disordered concurrence, we also discuss the 
results for the disordered single and two site entanglement entropy since by calculating 
the relevant disordered correlation functions for the concurrence, we already have all 
the necessary ingredients for the calculation. 

2. Entanglement Measures 

A number of entanglement measures exist, and each has advantages and disadvantages 
compared to the others. In this section, we briefly discuss two different entanglement 
measures to demonstrate how perturbation theory could be used to consider disorder 
in each of them. 

The concurrence, [291 [30] quantifies entanglement between two mixed qubits, 

C(p) = max{0, Ai - A 2 - A 3 - A 4 }, (2.1) 

where Ai > A 2 > A3 > A4. The Aj's are the square roots of the eigenvalues of the 
matrix p(a v (g) a v )p*(a v ® a v ), where a y is a Pauli spin matrix and p* is the complex 
conjugate of p. Concurrence has an advantage over other entanglement measures since 
is is relatively easy to calculate, and can be used for mixed states. Thus it can be 
used to study thermal entanglement. However, concurrence is limited in that it can 
only be used between pairs of qubits. 

The entropy of entanglement |31j measures the amount of entanglement in a pure 
state, pab- It is defined as the von Neumann entropy of a reduced density matrix, 
S (pa) = — tr (pa log pa) where pa = ^b(pab)- The entanglement entropy measures 
how mixed the subsystems, A and B, of a bipartite system are. Although this is 
only a measure of entanglement for pure states since for mixed states it measures 
both quantum and classical correlations, it has the advantage that it can be used to 
calculate entanglement between any two parts of a system. 

Since both of these entanglement measures use the density matrix p which can be 
written in terms of correlation functions, we can use many-body perturbation theory 
to find a perturbative expansion of each. For the concurrence, we only need two qubits, 
/ and m, of a system: 

Pl, m = £ (2.2) 

a,/3— 1,3;, y,z 

where (• • •) = tr(p---) is valid both for a thermal and a pure p. For two-site 
entanglement entropy, we need m = I + 1. The density matrix is even simpler if 
we wish to consider single-site entanglement entropy, then pa = J3 a =i x y z <J f( (J 'i)- 
Thus there are fewer correlation function averages to calculate. 

Other entanglement measures such as the relative entropy of entanglement, 
S(p\ \u css ) = — tr(plogp) — tr(plogcr css ) where o~ css is the closest separable state to p, 
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would be more difficult to calculate. In this case, we would need to find a css once the 
average had been taken. 

In the remainder of this paper, we concentrate on calculating the disorder 
averaged correlation functions necessary for finding the concurrence. 

3. The Model 

We consider an example, the XX spin chain in the thermodynamic limit, N —> oo, 
which means we can safely ignore boundary effects. This has Hamiltonian 

T N-l N 

i=i i=i 

where J is the coupling strength between neighbouring spins and B is an external 
magnetic field. Hq is the unperturbed part of the total Hamiltonian, H = Hq + Hi{ni) 
where Hi(fii) is the perturbation, and \ii is a random variable. We could use any 
distribution for the random term, but we consider the case when \xi = ji is taken from 
a Gaussian distribution centred at zero with variance A and 

1 N ^ 

Hi(ji) = -- 31 (of of +1 + ofof +1 ) . (3.2) 
i=i 

We apply a Jordan- Wigner transformation, ai = Ylm=i a rn ® (of + *of)/2 to the total 
Hamiltonian to get 

N-l n 
H = -J2(J + 3i) (aUl+l + a l+i a i) ~ B i 1 ' 2a ' a • ( 3 - 3 ) 

A Fourier transform, ai — J" ^e lql d(q) diagonalises the unperturbed Hamiltonian, 

leaving H = J* ^e(q)d^ (q)d(q) where e(q) = 2(B — Jcosq). Although we applied 
the Jordan- Wigner transformation to Hi , the Fourier transform is not useful. In order 
to treat Hi, we turn to fermionic many-body perturbation theory using equation [231 
The diagonal form of Hq allows us to calculate useful quantities such as the 
partition function, Zq = tr[exp(— /3Hq)], 

In Z = N f — In 1 2 cosh 
J-* 2tt \ 

where is the temperature and we have set ks = 1. 

4. Concurrence 

As discussed previously, to take disorder into account when calculating concur- 
rence, we must find the disorder average of each correlation function. For 
the XX spin chain many of the correlation functions are zero. The concur- 
rence between sites / and / + R in terms of correla tion functions is C = 
imax jo, \(<jfo-f +R + of <rf +R >| - ^(1 + (ofof +R )) 2 - 4(crf) 2 |. Thus for the disorder 

averaged concurrence, we need to calculate (crfaf +R + o-faf +R ), (o~?o-f +R ) and (of). 
The bar indicates the average over the disorder while the brackets, (• • •) = ti(p- ■ ■) 
denotes the thermal average. Each of these expectation values can be calculated 



(3.4) 
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as a perturbative expansion in A, the variance of the random distribution, using a 
functional integral technique. In this paper, we consider nearest and next nearest 
neighbour concurrence, where R = 1 and 2 respectively, however, the method we 
discuss applies to any value of R. 

The expectation values for the unperturbed XX chain were found in [32] . We 
now outline how they are calculated for the unperturbed system before discussing 
how to find their perturbative expansion. Using the Jordan- Wigner transformation, 
the necessary correlation functions are (o~fcrf +R + o-faf +R ) — 2(a\ IIm=i+i(^ — 
2a] n a m )ai +R + a\ +R ]\^2 l+1 {l-2al n a m )ai) 1 (afaf +R ) = ((l-2a]ai){l-2aj +R ai +R )} 

and (erf) = (1 - 2a\ai). 

Wick's theorem allows us to express n-point fermionic (and bosonic) correlation 
functions in terms of two-point correlation functions when Hq is a non-interacting 
system. Thus to calculate each of the above, we need only one quantity, G R = 
— (ajai+ii — a i a l+ R )i remembering the commutation relation, {aj, a m } — 5i. m . 

Applying the Fourier transform we used to diagonalise the unperturbed 
Hamiltonian, we get 

Gr = f g cos(qR)(l - 2d\q)d{q)). (4.1) 

When R — 0, this is equal to the magnetisation of a single site, Mq/N = Tln^i) /N . 
Magnetisation is defined as M — jj In Z and 

^/^gtanWfe)]. (4.2) 
Thus with the extra cos(qR) term in equation 14. 11 we have 

° R = f S C ° S{qR) tanh [PC{q)] ■ (4 ' 3) 

Once we consider disorder, again due to Wick's theorem (since our H is non- 
interacting), we must find an equivalent G R . Since all the correlation functions can 
be written in terms of G R , it is sufficient to calculate {d^ (q)d(q)) . 

5. Functional Integrals 

We introduce the functional partition function (or partition functional for convenience) 
which we will use to calculate the expectation values for the disorder averaged 
concurrence: 



Z [ry, rj\ = tr 



exp 




ill(r)ai +a\r)i{T) 



(5.1) 



Here, H is the Hamiltonian (equation Q3.3P ) and 77, 77 are Grassmann variables which 



We recover the partition function itself with Z[0, 0]. 
is the generating functional which we will use to 
calculate the correlation functions via functional derivatives. At the end of the 



are anti-commuting numbers [53 
The term J2i Vi( T ) a i + a 3^( T ) 
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calculation, we will set rj = rj = to regain the correct answer. Rewriting the above 
equation in functional integral form, we have 



J 23(7, 7 ) exp | -S + £ dr £ % (r) 7 j (r) + 7 (r)»» (r)] | , (5.2) 

where the action is 5 = f£ drj^i [7;( t )^t7;(t) + i?(7, 7)], and we have 
written a/ and in terms of Grassmann variables 7; and 7i- ^(7>7) = 

~ SHi 1 + ii) (77i+i + 7j+i7i) Sill (! ~ 2 77z) where we have suppressed the 
r argument after each term. Next we take n replicas of the partition functional, and 
average over the disorder ji 

z^wv\= {2n i )N/2 |n [* ex p (~^ 2 )] n z ^%]| , (5.3) 

where Z [rj al rj a ] gives a subscript a to every term in equation (|5.2|) . For each I, we can 
perform the average which corresponds to a Hubbard-Stratonovich transformation. 
This removes ji from Z n [rj, rj] leaving A. The result of this is that we swap the 
non-translationally invariant j\ term for a coupling between different replicas: 

W%^= ^ )7 )^J-^ + jfV^yTh,,(T)+7,>)^ (5.4) 



where S = S$ - AS? nt : 



S = dT Y] [%a 9 Tll,a ~ J (7,a7/+l,a + 7l + l,j!.a) + 2B7i,a7/,o] 

J ° La 

AS" t = y E {f dT ^2 [7i,,(T)7l+ 1 ,.(r) + 7i+i,a(T)7l,a(r)]| • (5.5) 

We can now diagonalise the unperturbed part of the Hamiltonian using a Fourier 
transform, 7J (t) = J_ dq/{2ir) exp(zgZ)7 a ("r, q). Applying this transformation to 
equation 15.41 we find 

= [ dr T ^-J2^(r,q)[dr+e(q)h a (r,q) 
Jo J-* 2vr ^ 

AS« t = £ T drdr' f f[ fe) J] 2*8{ qi q 2 + q 3 q,) (e<* + e~^) 

Z J ° J -* m =l V Z7T / o,a' 

x (e < «* + e- < «)7 (r 1 gi)7a(r,ga)7„'(r' l g3)7a'(^w)» ( 5 - 6 ) 

and dg/ (2tt) £ a [^ a ( T ) ?)7o(t, <?) + 7 a ( T , <?)?7a(T, q)], remembering that e{q) = 

2{B - Jcosq). 

In a system with no disorder and no source terms, we could now calculate the 
partition function, Zo[0, 0]. This calculation can either be completed in the path 
integral formalism or directly from the diagonalised Hamiltonian as found in equation 
13.41 We refer to, for example, pQ for the path integral calculation. 
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The next step is to define S' — \*f(d T + e)-f — rj-y — 777], where we have suppressed 
the equation's dependence on r, q and a for simplicity. Since 2? (7, 7) is invariant with 
respect to a translation, we let 7 — > 7 + / and 7 —¥ 7 + 5 where / and 5 are fields 
that we can choose the value of later. Substituting these new identities into S' and 
then expanding, we find that when (d T + e)f — r\ and g(— d T + e) = rj, many terms 
disappear, and we are left with 5" = 7 [d T + e(q)] 7 — rjf. Thus we need only find 
/ which we achieve by solving the inhomogeneous differential equation above. Using 
Green functions, /(r) = f£ dr'Q(T — t' , <?)?7(t') where 

G{t - r', q) = cxp [-e{q) (r - /)] {[1 - k(q)} 0(t - r') - k(q)e(r' - r)} , (5.7) 

and k(q) = (1 + e^ e ^)~ x . When r = t', we define the Green function as 
g(0-,q) = -k{q). Thus S' becomes dr £ £ a 7 (r, q) [d T + e(q)} 7o (r, 9) - 

/,f drc/V J*^ |f ^ a ^ (r, g)^(r - r', q)r] b (T' , q). The Green function term no longer 
depends on 7 and so can be taken out of the integral: 



Z n [i], 77] = cxp 



drdr' 



a, 6 



Z"[0,0]. 



(5.8) 



Where Z"[0,0] = / 2? (7, 7) exp + AS^ t ). We can also replace the 7s in the 
interaction term by taking advantage of functional derivatives: 



a^ 4 ^ f n (% ) e 2 * - * + «» - «*) ( gi 



m— 1 



+ e" 



2tt 



-2tt 



2tt- 



-2tt- 



^Vd( T ''1i) °~Vc(T',q 3 ) Si] b (T,q 2 ) 5r] a {T,qiy 



(5.9) 



Perturbation theory requires an expansion in a small term. For us, this is small 
term is A. Since A = jf, for the variance to be small, the average of the square 
of the random part of the coupling strength must also be small. Thus we expand 
exp(AS™ t ) = (1 + AS?, + (A^„ t ) 2 /2 + •••). We are left with 



Z" [rj, rj\ = cxp (-n In Z ) exp (ASSri) 

f drdr' [ p- y2rj a (T,q)g ab (T-T',q)r) b (T',q) 
Jo J-* 2tt ^ 



x cxp 



(5.10) 



Thus we have found a useful form of the disorder averaged replicated partition 
functional and are now in a position to calculate the correlation functions we require. 

The correlation functions are calculated by taking functional derivatives of the 
partition functional, and then setting rj = 77 = 0. 



(Tip(T 1 ,q 1 )ip(T 2 ,q2)) = lim 2ir- 



-2vr- 



n->0 5r ll (T 2 ,q 2 ) <frh( T i,gi) 
gives us a quenched disorder average, and 

1 - S S 



Z n [rj,v] 



rj—rj—0 



(Tip(T 1 ,q 1 )ip(T 2 ,q2)) 



-,2ir 



2tt- 



Z[rj,rj\ 5rn(T 2 ,q 2 ) ^(n.gi 



:Z[rj,v] 



(5.11) 



(5.12) 



rj—rj—O 
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Figure 1. The calculated Feynman diagrams: (a) Is the unperturbed diagram, 
and (b) and (c) are first order linked diagrams. Diagram (c) contributes only to 
the annealed average. 

the annealed disorder average. Here, T is the time ordering operator. Note that we 
need the average over the replicated partition functional, Z n \rj 1 rj\ for both, and the 
limit n — > in the quenched case means we are effectively calculating the average of 
the log of the partition functional. This is a variation of the replica method. We set 
n = 1 for the annealed case. 

In order to recover the needed correlation function, (dt (q)d(q)) , and recalling that 

WW, OW, Q')) = W, q)l(T',q>))9(T - t 1 ) - (7(r', q')j(T, <z)>0(r' - r), (5.13) 

we set t = t' and q — q' . When the rs are equal, the time ordering gives us 
— (7(1", q)"f{r, q)) which can be rewritten as —(d^(q)d(q)}. 

6. Perturbation Theory And The Linked Cluster Theorem 

Feynman diagrams |34j are an extremely useful tool which aid in the calculation 
of complicated functional derivatives by converting them into a diagrammatic form. 
These diagrams can be used to calculate the perturbative expansion of the correlation 
functions we require from Z n \fj, rj\ . When constructing the possible Feynman diagrams 
of the disorder averaged correlation functions, we will discover that both linked and 
vacuum diagrams exist. Vacuum diagrams are those which contain two or more parts, 
while a linked diagram is completely connected. Examples of linked diagrams are 
shown in figures QJb) an d ( c )- 

The linked cluster theorem demonstrates that when calculating correlation 
functions, vacuum diagrams cancel out and only the linked diagrams contribute. The 
numerator of the correlation functions is formed of a multiple of a sum of vacuum 
diagrams and a sum of linked diagrams while the denominator is a sum of vacuum 
diagrams only. The sums of vacuum diagrams therefore cancel, and we can write the 
correlation function as a sum over the linked diagrams. 

The linked cluster theorem also shows that the replica method is exact in 
perturbation theory [1]. Taking the quenched limit is equivalent to eliminating 
diagrams containing fermion loops (as shown in figure [T] (c)). Fermion loops arise 
from the average introducing an interaction between different replicas. Prior to the 
average, a diagram with a loop is actually a vacuum diagram. In the quenched limit, 
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these diagrams are no longer in the sum of linked diagrams, instead appearing in 
the sum of vacuum diagrams. Thus the quenched average leaves us a sum of linked 
diagrams without loops. 

Since we will calculate both the quenched and annealed averages, we must find 
and calculate all linked diagrams. 



7. Green Functions 



Equation 14.11 shows that the correlation function we need for the average is (cft(q)d(q)) 
which is found from equations 15.111 and 15.121 for the quenched an annealed averages 
respectively. To first order perturbation theory, the Feynman diagrams contributing 
to the correlation function are shown in figure [T] Using the standard rules for finite 
temperature many-body perturbation theory, we can calculate the contribution of each 
diagram (see, for example, [Ml 155] ). 

For example, diagram figure [T] (b) gives us 



2 f 



f drdr' ^ TT (^) V 5 ab 8 cd 2n8 ( qi - q 2 + q 3 - q 4 ) (e^ + e~^) 



drdr' 

/Q J —K m =l x """ ' abed 

x (e^ + e~^)2n5 (q fi - q 3 ) Q (t - r', q ) 5p c 2n6 (q 4 - Ql ) 

xQ(t' - r,q 4 ) Sad^nS (q 2 - q a ) G (t - T ai q 2 ) S ba , (7.1) 

where Q (t — r',q) was defined in equation 15.71 The multiple of 2 at the beginning 
indicates there are two topologically distinct ways in which we can make this diagram. 
In order to recover the appropriate correlation function, we set qp — q a and rp = r Q , 
and perform the analytically possible integrals. This leaves us with the cos 2 [(p + q)/2] 
term in Gr below. 

We find that the expansion of the correlation function gives 



Gr 



* dq 



cos (qR) { [1 - 2k(q)] 



p + q\ J k(q)-k(p) Pk(q)[l - k(q)] 



Hp) - £(<?)] 



e(p) - e{q) 



(7.2) 



/3 n cospcosqk(p) [1 — k(q)] k(q) 



As discussed previously, for the quenched average, n —y so the final term in the 
above equation is zero while for the annealed average, n = 1. 

If we take the limit j3 — > oo, we gain the zero temperature equivalent: 



Gr = r\ 



(1 + e i7vR ) sin I — I - 2 sin (aR) 



nR 



r7r/2 dp r dq ( [1 + co spcos g] [cos(qR) - cos(pi?)] 



Hp)-e{q)Y 
[1 — cospcosg] [cos(gi?) — e l ^ R cos(pR)] \ 

Hp) -e{<i)] 2 



(7.3) 
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Figure 2. Quenched nearest neighbour concurrence around the quantum critical 
point. The solid lines are no disorder, the longer dashes A = 0.0005 and the short 
dashes A = 0.001. The blue lines are T = 0, the cyan lines T = 0.01, the purple 
lines T = 0.02 and the orange lines T = 0.03. 



where v(x) = 2(B + J cos a;) and 



cos 




l (B/J) 



for B < J 
for B > J. 



We note that the quenched and annealed averages are identical in the T = case. 

We are now in a position to calculate the concurrence as described in previous 
sections. In terms of G#, the nearest neighbour concurrence (R = 1) is 



Ci = <0,G 1 \ 



1 + G 2 -gI) -4Go 



and the next nearest neighbour concurrence (R = 2) is 



C 2 = <^ 0, 



G 2 S~i 



9 2 

" ^0 _ U 2 



(7.4) 



(7.5) 



Then using equations 7.1 (and l7.3l for zero temperature), we can find the quenched 
and annealed averages of the concurrence. 

We also plan to briefly discuss the single-site entanglement entropy, S(pi) = 
-{(1+Go) log 2 [(l+Go)/2] + (l-Go) log 2 [(l-G )/2]}/2 which measures entanglement 
at zero temperature between a single site and the rest of the spin chain. In addition 
we discuss the two-site entanglement entropy, S(pi t i+i) = — a log a — clogc — d\og[(b + 

G?)/4, 



-?>2, 



d)/(b - d)] - &log(6 2 - d 2 ) where a = [(1 + G ) 2 - G^j/4, b = (1 - G 



c = [(1 — Go) 2 — G^]/4 and d = — Gi/2, which again can be used at zero temperature, 
and measures the entanglement between the two sites and the rest of the chain. 
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Figure 3. Annealed nearest neighbour concurrence around the quantum critical 
point. The solid lines are no disorder, the longer dashes A = 0.0005 and the short 
dashes A = 0.001. The blue lines are T = 0, the cyan lines T = 0.01, the purple 
lines T = 0.02 and the orange lines T = 0.03. 

8. Disorder Averaged Concurrence 

We have restricted the plots to be close to the quantum critical point (QCP) which 
for this system occurs at B = J since this is the region where the effect of the disorder 
is the most interesting. 

Figure [2] demonstrates how the quenched nearest neighbour concurrence behaves 
as the magnetic field varies for different values of temperature and disorder. We find 
that at zero temperature, increasing A always decreases the concurrence, and no new 
entanglement is created above the QCP. For finite temperatures, there is a crossover 
point of the magnetic field, B c , below which increasing A decreases concurrence, 
and above which, increasing A increases concurrence. As the temperature increases, 
B c generally decreases while the amount that A increases the concurrence above B c 
decreases. 

For the annealed nearest neighbour concurrence, shown in figure [31 we find a 
similar pattern of behaviour. However, B c in this case is much lower than for 
the quenched average. Both the quenched and annealed next nearest neighbour 
concurrence, shown in figures |4] and [5] respectively again give similar results, but 
with lower values for the concurrence overall. 

A possible reason for the behaviour described above could be that the disorder 
has a similar effect on entanglement to the temperature. Both have the effect of 
mixing energy levels allowing for the possibility of creating entanglement above the 
QCP. However, at higher temperatures and lower magnetic field, the likely effect of 
more mixing is to decrease entanglement. 

A puzzling feature of these results is that disorder does not increase the 
concurrence above the QCP at zero temperature as we would expect following the 
reasoning in the argument above. We expect that we would need to calculate the 
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perturbation series to a higher order to see this behaviour. Alternatively, perhaps the 
value of A simply needs to be larger than perturbation theory is valid for to show this. 

We have also calculated the correlation functions for a random magnetic held, 
when Hi(bi) = ^2ibiaf. Again we find the same behaviour for the concurrence, 
with a crossover point, B c , close to the QCP for the quenched average and lower for 
the annealed average, and zero concurrence at T = for B > J for both nearest 
and next nearest neighbour concurrence. Close to the QCP, increasing A for bi, 
increases the concurrence more than for j% for both nearest and next nearest neighbour 
entanglement. 

In addition to the concurrence we have considered the single and two site 
entanglement entropy. Since Go = Go at zero temperature, disorder doesn't affect 
the single-site entanglement entropy in the regime we are able to apply perturbation 
theory in. This is the case for whether the randomness is in ji or bi . 

Interestingly, the two-site entanglement entropy actually increases with disorder 
for both ji and bi for all values of B. Increasing A increases the amount of 
entanglement by a larger amount the closer B gets to J for each. The random magnetic 
field, bi, has a greater effect on the entanglement that ji; increasing A by the same 
amount increases the entanglement entropy more for bi than for ji . 

We compare our results to [TT] which uses sampling to look at an XX spin 
chain in a random magnetic field. They consider concurrence at zero temperature, 
and the random field taken from a Gaussian distribution as well as a Lorentzian 
distribution. They find that increasing the disorder below the QCP decreases the 
nearest neighbour concurrence, while above the QCP, increasing disorder increases 
the concurrence. While our results for zero temperature agree with this behaviour 
below the QCP, above B = J, concurrence remains zero for us. However, we consider 
weak disorder while the lowest disorder [11] calculates is an order of magnitude larger 
than ours. It is possible that calculating the perturbative expansion to higher order 
would allow us to observe this behaviour. 

9. Conclusions 

We have found that for weak disorder, concurrence, the entanglement between 
qubit pairs, in general decreases with increasing quenched disorder. At non-zero 
temperatures, close to the QCP, disorder instead acts to increase entanglement. In the 
case of two-site entanglement entropy, (remembering this is a measure of entanglement 
only at zero temperature), the entanglement between two sites and the remainder of 
the spin chain, increasing the disorder increases entanglement. 

Our results demonstrate that disorder is not necessarily detrimental to 
entanglement and thus to schemes which use it as a resource. This is, however, 
dependent on which entanglement measure is appropriate and the parameter regime 
we consider. 

We can use the perturbation theory method to find the disorder averaged 
concurrence for other systems. However, if the system cannot be expressed as a free 
fermion (non-interacting) model, the number of disorder averaged correlation functions 
required would be increased since Wick's theorem could no longer be applied. For 

2 2 

example, we would no longer have (crfo-^ +1 ) = G — G Rl and would instead need to 

calculate a four-point correlation function since (ofaf +l ) = 2Go — \ + A(a\aia\ +R ai+R) . 
One strength of our paper is that we consider finite temperature as well as 



Concurrence in Disordered Systems 



13 




0.96 0.98 1.00 1.02 



Figure 4. Quenched next nearest neighbour concurrence around the quantum 
critical point. The solid lines are no disorder, the longer dashes A = 0.0005 and 
the short dashes A = 0.001. The blue lines are T = 0, the cyan lines T = 0.01, 
the purple lines T = 0.02 and the orange lines T = 0.03. 
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Figure 5. Annealed next nearest neighbour concurrence around the quantum 
critical point. The solid lines are no disorder, the longer dashes A = 0.0005 and 
the short dashes A = 0.001. The blue lines are T = 0, the cyan lines T = 0.01, 
the purple lines T = 0.02 and the orange lines T = 0.03. 

zero temperature entanglement. Since we use perturbative methods, our work has 
the propensity to be extended in many directions. For example, higher order 
calculations and application of the Dyson equation would be useful. In particular, 
perturbation theory allows for the consideration of time dependent non-equilibrium 
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disorder averaged measures of entanglement using techniques such as the Keldysh 
formalism. This would be interesting since it would show how disorder affects the 
finite temperature entanglement of a system as it evolves over time. 



References 



Altland A and Simons B 2006 Condensed matter field theory (Cambridge University Press) 
Abrikosov A A, Gor'kov L P and Dzyaloshinskii I. Ye. 1992 Methods of Quantum Field Theory 

in Statistical Physics (Pergamon Press) 
De Dominicis C and Giardina I 2006 Random Fields and Spin Glasses: A Field Theory Approach 

(Cambridge University Press) 
Dasgupta C and Ma S 1980 Phys. Rev. B 22 1305 
Fisher D S 1994 Phys. Rev. B 50 3799 

Amico L, Fazio R, Osterloh A and Vedral V 2008, Rev. Mod. Phys. 80 517 
Refael G and Moore J E 2004 Phys. Rev. Lett. 93, 260602 
Laflorencie N 2005 Phys. Rev. B 72, 140408(R) 
Santachiara R 2006 J. Stat. Mech. L06002 

De Chiara G, Montangcro S, Calabrese P and Fazio R 2006 J. Stat. Mech. P03001 
Fujinaga M and Hatano N 2007 J. Phys. Soc. Jpn. 76, 094001 

Binosi D, De Chiara G, Montangero S and Recati A 2007 Phys. Rev. B, 76, 140405(R) 
Refael G and Moore J E 2007 Phys. Rev. B 76, 024419 
Hoyos J A and Rigolin G 2006 Phys. Rev. A 74 062324 
Hide J, Son W and Vedral V 2009 Phys. Rev. Lett. 102 100503 
Hide J and Vedral V 2010 Physica E 42, 359 
Bose S 2003 Phys. Rev. Lett. 91 207901 
Bosc S 2007 Contemp. Phys. 48 13 

Bayat A and Karimipour V 2005 Phys. Rev. A 71 042330 
Bayat A and Bosc S 2010 Adv. Math. Phys. 2010 1 

Wang X, Bayat A, Schirmer S G and Bose S 2010 Phys. Rev. A 81 032312 
Burgarth D and Bose S 2005 New J. Phys. 7 135 
Burrell C K and Osborne T J 2007 Phys. Rev. Lett. 99 167201 
Burrell C K, Eisert J and Osborne T J 2009 Phys Rev A 80 052319 
Allcock J and Linden N 2009 Phys. Rev. Lett. 102 110501 

Petrosyan D, Nikolopoulos G M and Lambropoulos P 2010 Phys. Rev. A 81 042307 
Wells B O, Lee Y S, Kastner M A, Christianson R J, Birgcncau R J, Yamada K, Endoh Y and 

Shiranc G 1997 Science 277 1067 
Anicich P G O and Grinberg H 2002 Int. J. Quant. Chem. 90 1562 
Wootters W K 1998 Phys. Rev. Lett. 80 2245 
Hill S and Wootters W K 1997 Phys. Rev. Lett. 78 5022 

Bennett C H, Bernstein H J, Popescu S and Schumacher B 1996 Phys. Rev. A 53 2046 
Barouch E and McCoy B M 1971 Phys. Rev. A 3 786 
Le Bellac M 2000 Thermal Field Theory (Cambridge University Press) 

Mattuck R D 2006 A Guide to Feynman Diagrams in the Many-body Problem (Dover 
Publications) 

[35] Ncgclc J W and Orland H 1998 Quantum many-particle systems (Westview Press) 



